Anomalous sound detection training apparatus, acoustic feature extraction apparatus, anomalous sound sampling apparatus, and methods and programs for the same

ABSTRACT

An anomalous sound detection training apparatus includes: a first acoustic feature extraction unit that extracts an acoustic feature of normal sound based on training data for normal sound by using an acoustic feature extractor; a normal sound model updating unit that updates a normal sound model by using the acoustic feature extracted; a second acoustic feature extraction unit that extracts an acoustic feature of anomalous sound based on simulated anomalous sound and extracts the acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor; and an acoustic feature extractor updating unit that updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted, in which processing by the units is repeatedly performed.

TECHNICAL FIELD

This invention relates to techniques for detecting anomalous sound of a machine from sound signals. For instance, it relates to techniques for generating an acoustic feature extractor for anomalous sound detection.

BACKGROUND ART

In a factory or the like, even a shutdown of an industrial device installed therein, such as a large-sized manufacturing or molding machine, due to a failure significantly hinders the factory's operation. It is thus necessary to routinely monitor the operational status of the device and immediately take measures in the event of an anomaly. One solution is to regularly dispatch a maintenance person to the field from an industrial device management service to check machine components for wear and the like. Since this requires enormous personnel and/or traveling costs and labor, however, it is difficult to implement this on all industrial devices and/or factories.

A solution to this is to install a microphone inside a machine so as to routinely monitor the operational sound of the machine. Through analysis of the operational sound, any occurrence of sound that is likely to be an anomaly (i.e., anomalous sound) is detected and an alert is raised, thereby solving the anomaly. However, setting the types of anomalous sound and/or methods for detecting them for every machine type or individual unit is even more expensive than manual monitoring. Thus, there is a need for automated design of rules for automatically detecting anomalous sound.

As a way to address this problem, anomalous sound detection based on statistical approaches is well known (see Non-patent Literature 1, for instance). Anomalous sound detection based on statistical approaches is generally classified into supervised anomalous sound detection and unsupervised anomalous sound detection. In supervised anomalous sound detection, a classifier is trained from training data for normal sound and anomalous sound; whereas in unsupervised anomalous sound detection, a classifier is trained only from training data for normal sound. For an industrial application, unsupervised anomalous sound detection is often used because training data for anomalous sound is difficult to collect.

A training/detection flow for unsupervised anomalous sound detection is as shown in FIG. 9. During training, an acoustic feature obtained from sound data (training data) at the time of normal operation is extracted. Thereafter, a normal sound model (a probability density function) is trained from the acoustic feature. Then, during determination, an acoustic feature is extracted in relation to a newly obtained observation, and a negative logarithmic likelihood (i.e., degree of anomaly) is evaluated against an already trained normal sound model. If the value is smaller than a threshold, it is determined to be normal; and if the value is larger than the threshold, it is determined to be anomalous. In other words, this is evaluation of how well an observed sound fits the normal sound model. It is based on the idea that sound “resembling” training data for normal sound should be produced if the observation represents normal sound and sound “not resembling” the training data for normal sound should be produced if the observation represents anomalous.

A challenge in unsupervised anomaly detection is design of an acoustic feature extractor F(·). In supervised anomalous sound detection, an acoustic feature that can correctly identify a target of determination is manually designed. For example, if it is known that the normal sound is a sinusoidal wave of 1000 Hz and the anomalous sound is a sinusoidal wave of 2000 Hz, a log power of a Mel filter bank (log-MFBO) is extracted per frame because the two sounds have different tone timbres. If the normal sound is steady engine noise and the anomalous sound is a “bumping” sound emitted by hitting of devices, a temporal difference of the power of a Mel filter bank (ΔMFBO) is extracted because the anomalous sound is a sporadic sound. “Deep training”, a type of supervised training, is said to be able to automatically design an acoustic feature from training data.

PRIOR ART LITERATURE Non-Patent Literature

Non-patent Literature 1: Tsuyoshi Ide and Masashi Sugiyama, “Anomaly Detection and Change Detection”, Kodansha, pp. 6-7, 2015.

SUMMARY OF THE INVENTION Problems to be Solved by the Invention

In unsupervised anomaly detection, however, it cannot be known anomalous sound with what kind of sound characteristics will occur. It is accordingly difficult to manually design the acoustic feature extractor and the use of deep training is also difficult. For example, if the anomalous sound is assumed to be a sinusoidal wave of 2000 Hz and log-MFBO is used as the acoustic feature due to the fact that the normal sound was a sinusoidal wave of 1000 Hz, anomalous sound like a “bumping” sound emitted by hitting of devices cannot be detected. The opposite is also true. This has forced use of, for example, Mel filter bank cepstrum coefficient (MFCC), which is a general-purpose acoustic feature, resulting in low accuracy of detection compared to supervised training.

An object of the present invention is to provide an anomalous sound detection training apparatus that is capable of generating an acoustic feature extractor for anomalous sound detection irrespective of whether training data for anomalous sound is available or not, as well as an acoustic feature extraction apparatus, an anomalous sound sampling apparatus, methods, and programs associated with the anomalous sound detection training apparatus.

Means to Solve the Problems

An anomalous sound detection training apparatus according to an aspect of this invention includes: a first acoustic feature extraction unit that extracts an acoustic feature of normal sound based on training data for normal sound by using an acoustic feature extractor; a normal sound model updating unit that updates a normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit; a second acoustic feature extraction unit that extracts an acoustic feature of anomalous sound based on simulated anomalous sound and extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor; and an acoustic feature extractor updating unit that updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit. Processing by the first acoustic feature extraction unit, the normal sound model updating unit, the second acoustic feature extraction unit, and the acoustic feature extractor updating unit is repeatedly performed.

An anomalous sound detection training apparatus according to an aspect of this invention includes: a first acoustic feature extraction unit that extracts an acoustic feature of normal sound based on training data for normal sound by using an acoustic feature extractor; a normal sound model updating unit that updates a normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit; a second acoustic feature extraction unit that extracts an acoustic feature of anomalous sound based on training data for anomalous sound and extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor; and an acoustic feature extractor updating unit that updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit. Processing by the first acoustic feature extraction unit, the normal sound model updating unit, the second acoustic feature extraction unit, and the acoustic feature extractor updating unit is repeatedly performed.

Effects of the Invention

It is possible to generate an acoustic feature extractor for anomalous sound detection irrespective of whether training data for anomalous sound is available or not.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram for describing an example of an anomalous sound detection training apparatus.

FIG. 2 is a flowchart for describing an example of an anomalous sound detection training method.

FIG. 3 is a block diagram for describing an example of an anomalous sound detection apparatus.

FIG. 4 is a flowchart for describing an example of an anomalous sound detection method.

FIG. 5 is a diagram for describing an example of the structure of the present invention.

FIG. 6 is a diagram for describing an example of an acoustic feature function.

FIG. 7 shows an example of occurrence frequency distribution of an amplitude spectrum spectrum for a single sound signal.

FIG. 8 shows an example of sparsity of an amplitude spectrum.

FIG. 9 is a diagram for describing a conventional technique.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[Overview of the Invention]

[[Structure of the Invention]]

A hierarchical structure of the present invention is shown in FIG. 5.

(I) Irrespective of supervised/unsupervised, a “normal sound model” and an “acoustic feature extractor” are optimized simultaneously from training data.

(II) In the case of unsupervised, training data for anomalous sound is sampled.

(III) Irrespective of supervised/unsupervised, the “acoustic feature extractor” is represented as a matrix product and a logarithm operation, or as a multilayer perceptron, for example.

(IV) A training/detection algorithm that simultaneously implements (I) to (III) is provided.

The following will describe (I) to (III). Note that for the sake of description, they will be described in the order of (III), (I), (H). It will be followed by description on [Anomalous sound detection training apparatus and method] for performing them.

[[(III) Definition of Acoustic Feature Extractor]]

A series derived by Fourier-transformation of an observation signal and extraction of a power spectrum or amplitude spectrum is expressed as X_(ω,t)øR₊ ^(Ω×T). Here, ω={1, 2, . . . , Ω} and t={1, 2, . . . , T} represent indices of frequency and time, respectively, Ω represents the number of frequency bins, and T represents the number of observed frames. R₊ represents a set of non-negative real numbers.

Most cases of acoustic feature extraction can be represented as a matrix operation on this series with a logarithmic function. For example, log-MFBO can be described as spectrum X_(t)=(X_(1,t), . . . , X_(Ω,t))^(T) of all frequencies in a certain frame multiplied by M Mel filter banks mm, of which the log is then taken.

$\begin{matrix} {{F\left( {M,X_{t}} \right)} = {\begin{bmatrix} {\ln \left( {\max \left( {{m_{1}^{T}X_{t}},ɛ} \right)} \right)} \\ \vdots \\ {\ln \left( {\max \left( {{m_{M}^{T}X_{t}},ɛ} \right)} \right)} \end{bmatrix} = {\ln \left( {\max \left( {{MX}_{t},ɛ} \right)} \right)}}} & (1) \end{matrix}$

Here, T represents transposition. Also, when ΔMFBO is defined as the temporal difference of Mel filter bank's output, it can be described as follows using a difference matrix D.

$\begin{matrix} {{{F\left( {M,D,X_{t}} \right)} = {\ln \left( {\max \left( {{D\begin{bmatrix} {MX}_{t - 1} \\ {MX}_{t} \end{bmatrix}},ɛ} \right)} \right)}},{D = \begin{bmatrix} {- E_{M}} \\ E_{M} \end{bmatrix}^{T}}} & (2) \end{matrix}$

Here, E_(M) is a unit matrix of M dimensions.

Herein, an acoustic feature extractor is defined with extension of this concept. First, past P₁ frame through future P₂ frame are taken from determination target frame t, and a vertical connection of them is defined as x_(t), where x_(t) is represented as shown in FIG. 6 and the formula below.

$\begin{matrix} \begin{matrix} {x_{t} = \left( {\underset{\underset{X_{t - P_{1}}}{}}{X_{1,{t - P_{1}}},\ldots \mspace{14mu},X_{\Omega,{t - P_{1}}}},\underset{\underset{X_{{t - P_{1} + 1}\;}}{}}{X_{1,{t - P_{1} + 1}},\ldots \mspace{14mu},X_{\Omega,{t - P_{1} + 1}}},\ldots \mspace{14mu},\underset{\underset{X_{t + P_{2}}}{}}{X_{1,{t + P_{2}}},\ldots \mspace{14mu},X_{\Omega,{t + P_{2}}}}} \right)^{T}} \\ {= {\left( {x_{t,1},\ldots \mspace{14mu},x_{t,Q}} \right)^{T} \in R_{+}^{Q}}} \end{matrix} & \left( 1^{\prime} \right) \end{matrix}$

Here, Q=Ω×(P₁+P₂+1) holds. The modification to the first to second lines of formula (1′) is rewriting for the purpose of simplifying the notation; the contents and order of the elements are unchanged. An operation that multiplies x_(t) by a matrix A∈R^(D×Q) and takes the log thereof is defined as the acoustic feature extractor. That is, it will be as follows.

$\begin{matrix} {{F\left( {A,x_{i}} \right)} = {\begin{bmatrix} {f\left( {a_{1},x_{i}} \right)} \\ \vdots \\ {f\left( {a_{D},x_{i}} \right)} \end{bmatrix} = {\begin{bmatrix} {\ln \left( {\max \left( {{a_{1}^{T}x_{i}},ɛ} \right)} \right)} \\ \vdots \\ {\ln \left( {\max \left( {{a_{D}^{T}x_{i}},ɛ} \right)} \right)} \end{bmatrix} = {\ln \left( {\max \left( {{Ax}_{i},ɛ} \right)} \right)}}}} & (3) \\ {\mspace{79mu} {A = {\begin{bmatrix} a_{1,1} & \ldots & a_{1,Q} \\ \vdots & \ddots & \vdots \\ a_{D,1} & \ldots & a_{D,Q} \end{bmatrix} = \begin{bmatrix} a_{1} \\ \vdots \\ a_{D} \end{bmatrix}}}} & (4) \end{matrix}$

Thus, optimization of the acoustic feature extractor can be said to be optimization of matrix A.

[[(I) Simultaneous Optimization of Normal Sound Model and Acoustic Feature Extractor]]

First, a normal sound model (a probability density function followed by a feature amount extracted from normal sound) is defined as p(F(A,x)|z=0). Then, a negative logarithmic likelihood (the degree of anomaly) can be computed as below.

g(A,x)=−ln p(F(A,x)|z=0)  (5)

Here, although the normal sound model may be any of various types, when it is a Gaussian mixture distribution, for example, it can be described as:

$\begin{matrix} {{p\left( {{{F\left( {A,x} \right)}z} = 0} \right)} = {\sum\limits_{k = 1}^{K}{W_{k}\exp \left\{ {{- \frac{1}{2}}\left( {{F\left( {A,x} \right)} - \mu_{k}} \right)^{T}{\sum\limits_{k}^{- 1}\left( {{F\left( {A,x} \right)} - \mu_{k}} \right)}} \right\}}}} & (6) \\ {\mspace{79mu} {W_{k} = \frac{w_{k}}{\left( {2\pi} \right)^{D/2}\sqrt{\det \left( \sum\limits_{k} \right)}}}} & (7) \end{matrix}$

Here, K is the number of mixtures, and k, Σ_(k), and w_(k) are a mean vector, a covariance matrix, and a weight parameter for the kth distribution, respectively.

In anomalous sound detection, when the value of g(A,x) is greater than threshold ϕ, it is determined to be “anomalous”; when the value is smaller than threshold ϕ, it is determined to be “normal”. Although this determination can be described as a step function, for the simplicity of the subsequent calculations, it is described with a sigmoid function:

$\begin{matrix} {{\sigma \left\{ \alpha \right\}} = \frac{1}{1 + {\exp \left( {- \alpha} \right)}}} & (8) \end{matrix}$

The sigmoid function 0≤σ{α}≤1 outputs a value closer to 0 as α is smaller and outputs a value closer to 1 as α is greater. That is, when σ{g(A,x-ϕ)} is equal to or greater than 0.5, it is determined to be anomalous; when it is equal to or smaller than 0.5, it is determined to be normal.

An aim of an anomalous sound detection algorithm is to output a great a {g(A,x−ϕ)} when the observation represents anomalous sound and to output a small σ{g(A,x−ϕ)} when it represents normal sound. That is, minimizing the following formula is the aim of optimization of anomalous sound detection.

L=∫σ{g(A,x)−ϕ}(p(F(A,x)|z=1)−p(F(A,x)|z=0))dx  (9)

Here, p(F(A,x)|z=1) is an anomalous sound model (a probability density function followed by a feature amount extracted from anomalous sound).

From the foregoing, simultaneous optimization of the normal sound model and the acoustic feature extractor should be performed so as to simultaneously satisfy both:

$\begin{matrix} {{p\left( {{{F\left( {A,x} \right)}z} = 0} \right)} = {\underset{p{({{{F{({A,x})}}z} = 0})}}{\arg \; \min}L}} & (10) \\ {A = {\underset{A}{\arg \; \min}\mspace{14mu} L}} & (11) \end{matrix}$

Optimization of the acoustic feature extraction can be said to be design of the matrix A such that formula (11) is optimized.

Optimization of formula (10) can be carried out with a maximum likelihood method, for example. Thus, consider the optimization of formula (11). The foremost challenge in calculating formula (9) is the integration operation. Thus, this integration is replaced with arithmetic mean of training data. First, an objective function is defined as:

L=−∫σ{g(A,x)−ϕ}p(F(A,x)|z=1)dx+∫σ{g(A,x)−ϕ}p(F(A,x)|z=0)dx  (12)

Each of the terms can then be regarded as an expectation operation for σ{G(A,x−ϕ)}, thus they can be described as:

L=−E[{g(A,x)−ϕ}]_(x|z=1) +E[σ{g(A,x)−ϕ}]_(x|z=0)  (13)

Here, if E[·]_(x|z=0) is an expectation operation relating to normal sound, while E[·]_(x|z=1) is an expectation operation relating to anomalous sound.

Here, if replacing the expectation operations with arithmetic means, formula (11) can be approximated as follows:

$\begin{matrix} {A = {\underset{A}{\arg \; \min} - {\frac{1}{N}{\sum\limits_{n = 1}^{N}{\sigma \left\{ {{g\left( {A,x_{n}} \right)} - \varphi} \right\}}}} + {\frac{1}{T}{\sum\limits_{t = 1}^{T}{\sigma \left\{ {{g\left( {A,x_{t}} \right)} - \varphi} \right\} \mspace{50mu} (14)}}}}} \\ {= {\underset{A}{\arg \; \min}\underset{\underset{J}{}}{{{- \frac{1}{N}}{\sum\limits_{n = 1}^{N}\Psi_{n}}} + {\frac{1}{T}{\sum\limits_{t = 1}^{T}\Psi_{t}}}}\mspace{335mu} (15)}} \end{matrix}$

Here, N is a time index of training data for anomalous sound. That is, the accuracy of anomalous sound detection will be maximized by calculating the arithmetic means of Ψ_(1, . . . , T) using training data for normal sound, calculating the arithmetic means of Ψ_(1, . . . , T) using training data for anomalous sound, and designing A so that the sum of the two arithmetic means is minimized. In the case of supervised anomalous sound detection where training data for anomalous sound has been obtained, this can be performed without modification. In the case of unsupervised anomalous sound detection, N pieces of anomalous sound data are sampled (simulated) in accordance with the sampling algorithm described below. In that case, N should be set about N=T.

Now that the objective function has been set with formula (15), it is then minimized. For the minimization, a gradient method may be used. When a steepest descent method is used as the gradient method, the matrix may be optimized for each dimension d as:

$\begin{matrix} \left. a_{d}\leftarrow{a_{d} - {\lambda \frac{\partial J}{\partial a_{d}}}} \right. & (16) \end{matrix}$

Here, λ is step size. For this optimization, AdaGrad, that is, an improved algorithm of the steepest descent method, or the like may be used, of course.

In a case where the normal sound model is represented as a Gaussian mixture distribution, a gradient vector ∂J/∂a_(d) can be calculated follows using the chain rule of partial differentiation.

$\begin{matrix} {\frac{\partial J}{\partial a_{d}} = {{{- \frac{1}{N}}{\sum\limits_{n = 1}^{N}{\frac{\partial\Psi_{n}}{\partial{g\left( {A,x_{n}} \right)}}\frac{\partial{g\left( {A,x_{n}} \right)}}{\partial{F\left( {A,x_{n}} \right)}}\frac{\partial{F\left( {A,x_{n}} \right)}}{\partial a_{d}}}}} + {\frac{1}{T}{\sum\limits_{t = 1}^{T}{\frac{\partial\Psi_{t}}{\partial{g\left( {A,x_{t}} \right)}}\frac{\partial{g\left( {A,x_{t}} \right)}}{\partial{F\left( {A,x_{t}} \right)}}\frac{\partial{F\left( {A,x_{t}} \right)}}{\partial a_{d}}}}}}} & (17) \end{matrix}$

Here, each of the components can be calculated by:

$\begin{matrix} {\mspace{79mu} {\frac{\partial\Psi_{i}}{\partial{g\left( {A,x_{i}} \right)}} = {\sigma \left\{ {{g\left( {A,x_{i}} \right)} - \varphi} \right\} \left( {1 - {\sigma \left\{ {{g\left( {A,x_{i}} \right)} - \varphi} \right\}}} \right)}}} & (18) \\ {\mspace{79mu} {\frac{\partial{g\left( {A,x_{i}} \right)}}{\partial{F\left( {A,x_{i}} \right)}} = {\sum\limits_{k = 1}^{K}\; {\Gamma_{k,i}\left\{ {\left( {{F\left( {A,x_{i}} \right)} - \mu_{k}} \right)^{T}\Sigma_{k}^{- 1}} \right\}}}}} & (19) \\ {\Gamma_{k,i} = \frac{W_{k}\exp \left\{ {{- \frac{1}{2}}\left( {{F\left( {A,x_{i}} \right)} - \mu_{k}} \right)^{T}{\Sigma_{k}^{- 1}\left( {{F\left( {A,x_{i}} \right)} - \mu_{k}} \right)}} \right\}}{\sum\limits_{j = 1}^{K}\; {W_{j}\exp \left\{ {{- \frac{1}{2}}\left( {{F\left( {A,x_{i}} \right)} - \mu_{j}} \right)^{T}{\Sigma_{j}^{- 1}\left( {{F\left( {A,x_{i}} \right)} - \mu_{j}} \right)}} \right\}}}} & (20) \\ {\mspace{79mu} {\frac{\partial{F\left( {A,x_{i}} \right)}}{\partial a_{d}} = \begin{bmatrix} \frac{\partial{f\left( {a_{1},x_{i}} \right)}}{\partial a_{d,1}} & \ldots & \frac{\partial{f\left( {a_{1},x_{i}} \right)}}{\partial a_{1,Q}} \\ \vdots & \ddots & \vdots \\ \frac{\partial{f\left( {a_{D},x_{i}} \right)}}{\partial a_{d,1}} & \ldots & \frac{\partial{f\left( {a_{D},x_{i}} \right)}}{\partial a_{1,Q}} \end{bmatrix}}} & (21) \\ {\mspace{79mu} {\frac{\partial{f\left( {a_{j},x_{i}} \right)}}{\partial a_{d,q}} = \frac{C_{i,j,d,q}}{{ma}{\overset{.}{x}\left( {{a_{j}^{T}x_{i}},ɛ} \right)}}}} & (22) \\ {\mspace{79mu} {C_{i,j,d,q} = \left\{ \begin{matrix} {x_{i,q}{\delta \left( {d,j} \right)}} & \left( {{a_{j}^{T}x_{i}} > ɛ} \right) \\ 0 & ({otherwise}) \end{matrix} \right.}} & (23) \\ {\mspace{79mu} {{\delta \left( {d,j} \right)} = \left\{ \begin{matrix} 1 & {d = j} \\ 0 & ({otherwise}) \end{matrix} \right.}} & (24) \end{matrix}$

[[(II) Sampling of Training Data for Anomalous Sound]]

For calculation of formula (15) in unsupervised anomalous sound detection, how training data for anomalous sound is sampled will be discussed.

The distribution used with anomalous sound model p(F(A,x)|z=1) is arbitrary. In unsupervised training, no previous knowledge on anomalous sound has been obtained. In order to detect every kind of anomalous sound, it is desirable to construct an anomalous sound model with as least assumptions as possible. That is, it is desirable to sample anomalous sound only from generally known findings about the physical characteristics or frequency of appearance of an amplitude spectrum and/or power spectrum.

In the following, an anomalous sound model is constructed with assumptions made only for the occurrence frequency, sparsity, and continuity in temporal frequency direction of amplitude values.

First, consider the magnitude of X_(ω,n). FIG. 7 shows the occurrence frequency distribution of the amplitude spectrum for a single sound signal. The horizontal axis is strength and the vertical axis is the frequency of appearance. From FIG. 7, it can be seen that strengths around 0 are most likely to occur in the amplitude spectrum of a sound signal and the frequency of appearance exponentially decreases as the strength increases. X_(ω,n) should be sampled from a probability distribution that satisfies these characteristics.

A distribution that satisfies the characteristics is an exponential distribution, for example. An exponential distribution is a distribution such that:

$\begin{matrix} {{p\left( {o;\beta} \right)} = {\frac{1}{\beta}\exp \left\{ {- \frac{o}{\beta}} \right\}}} & (25) \end{matrix}$

In the distribution, the frequency of appearance of o decreases as the value of o increases. β is an average value. In the case of sampling X_(ω,n) from an exponential distribution, it will be as follows.

$\begin{matrix} {{\left. X_{\omega,n} \right.\sim\frac{1}{\beta}}\exp \left\{ {{- \frac{1}{\beta}}X_{\omega,n}} \right\}} & (26) \end{matrix}$

Since β is a parameter relating to the volume of observed sound, it can be set from the training data for normal sound. For example, it can be set as:

$\begin{matrix} {\beta = {\frac{1}{\Omega \; T}{\sum\limits_{\omega = 1}^{\Omega}\; {\sum\limits_{t = 1}^{T}\; X_{\omega,t}}}}} & (27) \end{matrix}$

If previous knowledge on stationary noise contained in observation signals or the like is available, this sampling may be modified as desired. For example, if an average spectrum Y_(ω) of stationary noise is known, sampling may be performed such that:

$\begin{matrix} {X_{\omega,n} = {Y_{\omega} + E_{\omega,n}}} & (28) \\ {{\left. E_{\omega,n} \right.\sim\frac{1}{\beta_{\omega}}}\exp \left\{ {{- \frac{1}{\beta_{\omega}}}E_{\omega,n}} \right\}} & (29) \end{matrix}$

β here can be calculated by:

$\begin{matrix} {\beta_{\omega} = \sqrt{\frac{1}{T}{\sum\limits_{t = 1}^{T}\; \left( {\max \left( {0,{X_{\omega,t} - Y_{\omega}}} \right)} \right)^{2}}}} & (30) \end{matrix}$

When sampled with formula (29), a volume average β_(ω) of anomalous sound will be constant independently of time. Depending on the distance to a microphone or the type of anomalous sound, the amplitude value, or the volume, might increase or decrease. Hence, it is desirable to vary the value of β_(ω,n) on a per-n basis. For representing this, formula (29) may be modified as follows, for example:

$\begin{matrix} {{\left. E_{\omega,n} \right.\sim\frac{1}{\beta_{\omega,n}}}\exp \left\{ {{- \frac{1}{\beta_{\omega,n}}}E_{\omega,n}} \right\}} & (31) \\ {\beta_{\omega,n} = {g_{\omega} \times b_{\omega}10^{\theta/20}}} & (32) \\ {g_{\omega} = \sqrt{\frac{1}{T}{\sum\limits_{t = 1}^{T}\; \left( {X_{\omega,t} - N_{\omega}} \right)^{2}}}} & (33) \\ {\left. b_{\omega} \right.\sim{{Uniform}\left( {0,1} \right)}} & (34) \end{matrix}$

Here, θ is a parameter in unit of dB for correcting the maximum volume of anomalous sound, and Uniform(a,b) is a continuous uniform distribution of maximum a and minimum b. For example, θ may be set to about 0 to 12. As another way, similar effects could be provided by sampling θ from a normal distribution or the like and setting it always to, for example, 1, without performing sampling of b_(ω).

It is also possible to focus on the sparsity of an amplitude spectrum (almost all of the elements being 0). FIG. 8 shows an example of plotting X_(1, . . . , Ω,t) for a certain frame. FIG. 8 indicates that a sound signal often has only one or several frequency components in a single frame. Focusing on this, the sampling of formula (31) may be modified as follows:

$\begin{matrix} {{\left. E_{\omega,n} \right.\sim\frac{1}{\beta_{\omega,n}}}\exp \left\{ {{- \frac{1}{\beta_{\omega,n}}}E_{\omega,n}} \right\}} & (35) \\ {\beta_{\omega,n} = \left\{ \begin{matrix} g_{\omega} & {w = z_{n}} \\ 0 & ({otherwise}) \end{matrix} \right.} & (36) \\ {\left. z_{n} \right.\sim{{Cat}\left( {1_{\Omega}/\Omega} \right)}} & (37) \end{matrix}$

Here, Cat(π) is category distribution and 1_(Ω) is a vector with a length of Ω whose elements are all 1. This makes only one spectrum bin become active at each time via z_(n), thereby driving the anomalous sound spectrum E to be sparse.

Also, sampling with formula (35) results in generation of a value of E_(ω,n) that is completely independent from the temporal frequency direction. A sound signal typically changes smoothly in terms of time, so that the spectrum thereof smoothly changes in the temporal frequency direction. For achieving this, F_(ω,n) may be smoothed in the temporal frequency direction.

While the smoothing can be done in various ways, it can be performed using a Gaussian filter as follows, for example:

X _(ω,n) =Y _(ω)+conv2(E _(ω,n) ,G(F,T))  (38)

Here, conv2 is a convolution operation of two dimensions, and G(F,T) is a two-dimensional Gaussian filter with a size of F×T. Although the parameters should be controlled experimentally, they may be set at about F=11 and T=31, for example.

[Anomalous Sound Detection Training Apparatus and Method]

As exemplarily shown in FIG. 1, an anomalous sound detection training apparatus includes a spectrum calculation unit 1, an initialization unit 2, a first acoustic feature extraction unit 3, a normal sound model updating unit 4, an anomalous sound sampling unit 5, a second acoustic feature extraction unit 6, an acoustic feature extractor updating unit 7, and a control unit 8. An anomalous sound detection training method is implemented by the units of the anomalous sound detection training apparatus executing the processing at step S1 to S82 described in FIG. 2 and below.

<Spectrum Calculation Unit 1>

Training data for normal sound is input to the spectrum calculation unit 1.

The training data for normal sound is sound signals collected by a microphone installed near (or alternatively inside) a device for which anomalous sound is to be detected. While the sampling frequency for collecting sound signals should be set so that the frequencies of normal sound can be recorded, it may be set at 16 kHz, for example.

The spectrum calculation unit 1 converts the training data into a frequency domain via an approach such as short-time Fourier transformation (STFT).

It then takes an absolute value of a complex spectrum to determine an amplitude spectrum X_(ω,t) (step S1). Alternatively, X_(ω,t) may be a power spectrum derived by determining the square root of the amplitude spectrum. Here, Fourier transformation length may be set at 32 ms and shift width may be set at 16 ms, for example.

In the following, X_(ω,t) as either an amplitude spectrum or a power spectrum will be abbreviated as spectrum X_(ω,t). The spectrum X_(ω,t) is output to the first acoustic feature extraction unit 3 and the second acoustic feature extraction unit 6.

If training data for anomalous sound is available, namely in the case of supervised training, the spectrum calculation unit 1 calculates the spectrum of anomalous sound based on the training data for anomalous sound. The spectrum of anomalous sound is the amplitude spectrum or the power spectrum of the anomalous sound, for example. The calculated spectrum of anomalous sound is output to the first acoustic feature extraction unit 3 and the second acoustic feature extraction unit 6.

<Initialization Unit 2>

Preset parameters for a training algorithm and preset parameters for an acoustic feature extractor are input to the initialization unit 2.

The various parameters may be set such that P₁=3, P₂=3, ε=min(X_(ω,t))/2, ϕ=100, and λ=5.0×10⁻³, for example.

The initialization unit 2 initializes the matrix A for the acoustic feature extractor and the normal sound model (step S2). After being initialized, the matrix A and the normal sound model are output to the first acoustic feature extraction unit 3.

The matrix A may be initialized with uniform random numbers, for example. While initialization of the normal sound model varies depending on the type of the probability density function, when a Gaussian mixture distribution is used, for example, the normal sound model can be initialized via a k-means algorithm or the like.

<First Acoustic Feature Extraction Unit 3>

To the first acoustic feature extraction unit 3, the acoustic feature function initialized by the initialization unit 2 and the training data for normal sound are input.

The first acoustic feature extraction unit 3 extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor (step S3). In this example, as the training data for normal sound, the first acoustic feature extraction unit 3 uses the spectrum X_(ω,t) calculated by the spectrum calculation unit 1 based on the training data for normal sound to extract the acoustic feature of normal sound. The extracted acoustic feature is output to the normal sound model updating unit 4.

A first round of processing by the first acoustic feature extraction unit 3 is performed by using the acoustic feature extractor initialized by the initialization unit 2. Second and subsequent rounds of processing by the first acoustic feature extraction unit 3 are performed by using the acoustic feature extractor as updated by the acoustic feature extractor updating unit 7.

If the spectrum of anomalous sound has been calculated by the spectrum calculation unit 1 based on training data for anomalous sound, the first acoustic feature extraction unit 3 extracts the acoustic feature of anomalous sound based on this calculated spectrum of anomalous sound by using the acoustic feature extractor.

<Normal Sound Model Updating Unit 4>

The normal sound model updating unit 4 updates the normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit 3 (step S4).

While the way of update varies depending on the type of the probability density function for the normal sound model, it may be updated by maximum likelihood estimation, for example. In a case where a Gaussian mixture distribution is used as the normal sound model, one to several steps of update may be performed in an EM algorithm (see Reference Literature 1, for instance).

-   [Reference Literature 1] Sadanori Konishi, “Introduction to     Multivariate Analysis, Appendix C: EM algorithm”, pp. 294-298,     Iwanami Shoten, 2010.

<Anomalous Sound Sampling Unit 5>

The anomalous sound sampling unit 5 samples simulated anomalous sound (step S5). The sampled anomalous sound is output to the second acoustic feature extraction unit 6.

The anomalous sound sampling unit 5 performs sampling of simulated anomalous sound based on formula (26), for example. The parameter β in formula (26) can be determined based on formula (27), for example. In this case, for X_(ω,t) in formula (27), the spectrum X_(ω,t) of normal sound calculated by the spectrum calculation unit 1 may be used.

If training data for anomalous sound is available, namely in the case of supervised training, sampling is not performed. That is, the subsequent processing may be performed using the training data for anomalous sound as the sampling result.

<Second Acoustic Feature Extraction Unit 6>

The second acoustic feature extraction unit 6 extracts the acoustic feature of anomalous sound based on simulated anomalous sound and extracts the acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor (step S6). The acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted are output to the acoustic feature extractor updating unit 7.

In this example, as the training data for normal sound, the second acoustic feature extraction unit 6 uses the spectrum X_(ω,t) calculated by the spectrum calculation unit 1 based on the training data for normal sound to extract the acoustic feature of normal sound.

The first round of processing by the second acoustic feature extraction unit 6 is performed by using the acoustic feature extractor initialized by the initialization unit 2. The second and subsequent rounds of processing by the second acoustic feature extraction unit 6 are performed by using the acoustic feature extractor as updated by the acoustic feature extractor updating unit 7.

<Acoustic Feature Extractor Updating Unit 7>

The acoustic feature extractor updating unit 7 updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit 6 (step S7). The updated acoustic feature extractor is output to the first acoustic feature extraction unit 3 and the second acoustic feature extraction unit 6. Also, the acoustic feature extractor as finally updated after repeated control by the control unit 8 is output as a final result of training performed by the anomalous sound detection training apparatus and method.

The acoustic feature extractor updating unit 7 updates the acoustic feature extractor based on formula (16), for example. Prior to execution of this update, threshold ϕ may be updated according to the formula below using training data for normal sound x₁, . . . , x_(T) and sampled anomalous sound data x₁, . . . x_(N). This stabilizes convergence.

$\begin{matrix} {\varphi = \frac{{{median}\left( {g\left( {A,x_{1,\ldots,T}} \right)} \right)} + {{median}\left( {g\left( {A,x_{1,\ldots,N}} \right)} \right)}}{2}} & (39) \end{matrix}$

<Control Unit 8>

The control unit 8 performs control so as to repeatedly perform processing by the first acoustic feature extraction unit 3, the normal sound model updating unit 4, the second acoustic feature extraction unit 6, and the acoustic feature extractor updating unit 7 (steps S81 and S82). This repeated processing is performed until the acoustic feature extractor and the normal sound model converge.

For example, the control unit 8 performs control so as to repeatedly perform the processing by the second acoustic feature extraction unit 6 and the acoustic feature extractor updating unit 7 (hereinafter referred to as a first repeated processing). Setting a first convergence determination condition as that the number of executions of processing has reached a certain number of times (for example, 50) or that an update norm for the matrix A has fallen below a certain value, the control unit 8 performs control so as to repeatedly perform the first repeated processing until the first convergence condition is satisfied, or in other words, until it is determined that convergence has been reached. The first convergence determination condition may also be a different condition.

The control unit 8 subsequently performs control so as to repeatedly perform the processing by the first acoustic feature extraction unit 3 and the normal sound model updating unit 4 and the first repeated processing (hereinafter referred to as a second repeated processing). Setting a second convergence condition as that the number of executions of processing has reached a certain number of times (for example, 200), the control unit 8 performs control so as to repeatedly perform the second repeated processing until the second convergence condition is satisfied, or in other words, until it is determined that convergence has been reached. The second convergence determination condition may also be a different condition.

[Anomalous Sound Detection Apparatus and Method]

As exemplarily shown in FIG. 3, an anomalous sound detection apparatus includes a spectrum calculation unit 11, an acoustic feature extraction unit 12, a degree-of-anomaly calculation unit 13, and a determination unit 14. An anomalous sound detection method is implemented by the units of the anomalous sound detection apparatus executing the processing at step S11 to S14 described in FIG. 4 and below.

<Spectrum Calculation Unit 11>

Operational sound of a machine for which anomalous sound is to be detected is collected through a microphone. The sampling rate used for this collection is similar to the one used for training. The collected sound signal is input to the spectrum calculation unit 11.

The spectrum calculation unit 11 obtains spectrum. X_(ω,t) based on the collected sound signal in a similar manner to the spectrum calculation unit 1 (step S11). The obtained spectrum X_(ω,t) is output to the acoustic feature extraction unit 12.

<Acoustic Feature Extraction Unit 12>

The acoustic feature extraction unit 12 uses the Obtained spectrum X_(ω,t) to extract the acoustic feature of the collected sound signal based on the acoustic feature extractor that has been output by the anomalous sound detection training apparatus and method as the final result of training (step S12). The extracted acoustic feature is output to the degree-of-anomaly calculation unit 13.

<Degree-of-Anomaly Calculation Unit 3>

The degree-of-anomaly calculation unit 13 calculates a negative logarithmic likelihood g(A,x) by using the extracted acoustic feature and a normal sound model that has been output by the anomalous sound detection training apparatus and method as the final result of training (step S13). Hereinafter, a negative logarithmic likelihood will be referred to as the degree of anomaly. The calculated degree of anomaly is output to the determination unit 14.

<Determination Unit 14>

The determination unit 14 outputs “anomalous” if the degree of anomaly of a current frame is equal to or greater than threshold ϕ (step S14). Although the threshold should be adjusted in accordance with the machine and/or environment in question, it may be set at about 1500, for example.

As with voice activity detection, “hangover”, which suppresses detection errors with heuristic rules, may be used. While the hangover processing to be applied may be any of various types, the hangover processing should be set in accordance with the types of false detection of anomalous sound.

As an example, musical noise that occurs during noise suppression could be determined as a sporadic anomalous sound. Since sound like a sporadic hitting sound often exhibits change in spectrum shape for 100 ms or more, the degree of anomaly would remain equal to or greater than the threshold for _(┌)(100/the frame shift width for STFT)_(┐) frames continuously. However, because in musical noise an anomalous amplitude spectrum value only occurs in the relevant frames, there would be several frames at most where the degree of anomaly remains equal to or greater than the threshold continuously. Thus, a rule for anomaly determination may be set as “outputting “anomalous” if the degree of anomaly remains equal to or greater than the threshold for F₁ frames or more continuously”, for example.

As another example, it is also conceivable that anomalous sound continues for a long period of time with the degree of anomaly slightly below the threshold due to a low volume of the anomalous sound. In such a case, a rule like “detecting a sound as anomalous if a total sum of its degree of anomaly over the last F₂ frames is equal to or greater than ϕ₁” may be added as a determination rule for continuous anomalous sound. Although ϕ₁ should be determined by tuning, it may be set at about ϕ₁=F₂×(ϕ−250), for example.

By detecting anomalous sound from a large-sized manufacturing or molding machine installed in a factory or the like using the anomalous sound detection apparatus and method as described above, faster handling of a failure and/or failure prediction become possible. This can contribute to increased efficiency of an industry, for example, manufacturing industry in particular.

[Program and Recording Medium]

In the case of implementing the processing in the anomalous sound detection training apparatus or the anomalous sound detection apparatus with a computer, the processing details of the functions to be provided by the anomalous sound detection training apparatus or the anomalous sound detection apparatus are described by a program. By the computer then executing the program, the processing is implemented on the computer.

The program describing the processing details may be recorded in a computer-readable recording medium. The computer-readable recording medium may be any kind of medium, for example, a magnetic recording device, an optical disk, a magneto-optical recording medium, and a semiconductor memory.

In addition, individual processing means may be embodied by execution of a predetermined program on a computer or at least some of their processing details may be implemented in hardware.

[Modifications]

Using the optimization metrics for acoustic feature described above, a multilayer perceptron can also be trained unsupervised. Specifically, although the acoustic feature extractor was described using formula (3), it may instead be implemented with a multilayer perceptron. For example, when a simplest fully-connected multilayer perceptron is used, the acoustic feature extractor would be as shown below. In this manner, the acoustic feature extractor may also be represented by a multilayer perceptron.

F(Xτ)=W ^((L)) g _(τ) ^((L-1)) +b ^((L))  (40)

g _(τ) ^((k))=σ_(θ) {W ^((k)) g _(τ) ^((k-1)) +b ^((k))}  (41)

Here, L, (k), and b(k) represent the number of layers of a neural network, a weight matrix, and a bias vector, respectively. L is a predetermined positive integer equal to or greater than 2. σ_(θ) is an activation function and may be implemented as a sigmoid function, for example. Input to the input layer may be such that g_(τ) ⁽¹⁾=x_(τ) or g_(τ) ⁽¹⁾=ln x_(τ), for example. g_(τ) ^((k)) is input to the kth layer of the multilayer perceptron. In this case, optimization of the acoustic feature extractor is optimization of Θ={(W(k),b(k)}, (k∈{2, . . . , L}).

For this optimization, a gradient method may be used as in formula (16). If a steepest descent method is used as the gradient method, it may be done such that:

$\begin{matrix} \left. W^{(k)}\leftarrow{W^{(k)} - {\lambda \frac{\partial J}{\partial W^{(k)}}}} \right. & (42) \\ \left. b^{(k)}\leftarrow{b^{(k)} - {\lambda \frac{\partial J}{\partial b^{(k)}}}} \right. & (43) \end{matrix}$

For this optimization, AdaGrad, that is, an improved algorithm of the steepest descent method, or the like may be used, of course. The objective function J is given by

$\begin{matrix} {J = {{{- \frac{1}{N}}{\sum\limits_{n = 1}^{N}\; \Psi_{n}}} + {\frac{1}{T}{\sum\limits_{t = 1}^{T}\; \Psi_{t}}}}} & (44) \end{matrix}$

Thus, formulae (42) and (43) can be described as:

$\begin{matrix} \left. W^{(k)}\leftarrow{W^{(k)} - {\lambda \left\lbrack {{{- \frac{1}{N}}{\sum\limits_{n = 1}^{N}\frac{\partial\Psi_{n}}{\partial W^{(k)}}}} + {\frac{1}{T}{\sum\limits_{t = 1}^{T}\frac{\partial\Psi_{t}}{\partial W^{(k)}}}}} \right\rbrack}} \right. & (45) \\ \left. b^{(k)}\leftarrow{b^{(k)} - {\lambda \left\lbrack {{{- \frac{1}{N}}{\sum\limits_{n = 1}^{N}\frac{\partial\Psi_{n}}{\partial b^{(k)}}}} + {\frac{1}{T}{\sum\limits_{t = 1}^{T}\frac{\partial\Psi_{t}}{\partial b^{(k)}}}}} \right\rbrack}} \right. & (46) \end{matrix}$

The gradients:

$\frac{\partial\Psi_{i}}{\partial W^{(k)}},\frac{\partial\Psi_{i}}{\partial b^{(k)}}$

can be calculated as follows using the chain rule of partial differentiation from formula (17).

$\begin{matrix} {\frac{\partial\Psi_{i}}{\partial W^{(k)}} = {\frac{\partial\Psi_{i}}{\partial{G\left( x_{i} \right)}}\frac{\partial{G\left( x_{i} \right)}}{\partial{F\left( x_{i} \right)}}\frac{\partial{F\left( x_{i} \right)}}{\partial W^{(k)}}}} & (47) \\ {\frac{\partial\Psi_{i}}{\partial b^{(k)}} = {\frac{\partial\Psi_{i}}{\partial{G\left( x_{i} \right)}}\frac{\partial{G\left( x_{i} \right)}}{\partial{F\left( x_{i} \right)}}\frac{\partial{F\left( x_{i} \right)}}{\partial b^{(k)}}}} & (48) \end{matrix}$

In a multilayer perceptron, gradients in the form of formulae (47) and (48) can be efficiently calculated by using error back propagation (see Reference Literature 2, for instance).

-   [Reference Literature 2] Hideki Asoh and seven others, “Deep     Learing”, the Japanese Society for Artificial Intelligence, Kindai     kagaku sha, pp. 5, 2015.

In addition to being executed chronologically in the order described, the processing described for the anomalous sound detection training apparatus or the anomalous sound detection apparatus may also be executed in parallel or individually depending on the processing ability of the apparatus executing the processing or on any necessity.

It goes without saying that other modifications may be made as desired without departing from the scope of this invention. 

1-14. (canceled)
 15. An anomalous sound detection training apparatus comprising: a first acoustic feature extraction unit that extracts an acoustic feature of normal sound based on training data for normal sound by using an acoustic feature extractor; a normal sound model updating unit that updates a normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit; a second acoustic feature extraction unit that extracts an acoustic feature of anomalous sound based on simulated anomalous sound and extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor; and an acoustic feature extractor updating unit that updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit, wherein processing by the first acoustic feature extraction unit, the normal sound model updating unit, the second acoustic feature extraction unit, and the acoustic feature extractor updating unit is repeatedly performed.
 16. An anomalous sound detection training apparatus comprising: a first acoustic feature extraction unit that extracts an acoustic feature of normal sound based on training data for normal sound by using an acoustic feature extractor; a normal sound model updating unit that updates a normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit; a second acoustic feature extraction unit that extracts an acoustic feature of anomalous sound based on training data for anomalous sound and extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor; and an acoustic feature extractor updating unit that updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit, wherein processing by the first acoustic feature extraction unit, the normal sound model updating unit, the second acoustic feature extraction unit, and the acoustic feature extractor updating unit is repeatedly performed.
 17. The anomalous sound detection training apparatus according to claim 15, further comprising: an anomalous sound sampling unit that samples the simulated anomalous sound, wherein the simulated anomalous sound is anomalous sound simulated by making assumption on one of occurrence frequency, sparsity, and continuity in temporal frequency direction of amplitude values.
 18. The anomalous sound detection training apparatus according to claim 17, wherein the simulated anomalous sound is X_(ω,n) as sampled based on a distribution: $\begin{matrix} {{\left. X_{\omega,n} \right.\sim\frac{1}{\beta}}\exp \left\{ {{- \frac{1}{\beta}}X_{\omega,n}} \right\}} & (26) \end{matrix}$ where ω is frequency, n is frame number, and β is an integer.
 19. An anomalous sound detection training method comprising: a first acoustic feature extraction step in which a first acoustic feature extraction unit extracts an acoustic feature of normal sound based on training data for normal sound by using an acoustic feature extractor; a normal sound model updating step in which a normal sound model updating unit updates a normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit; a second acoustic feature extraction step in which a second acoustic feature extraction unit extracts an acoustic feature of anomalous sound based on simulated anomalous sound and extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor; and an acoustic feature extractor updating step in which an acoustic feature extractor updating unit updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit, wherein processing at the first acoustic feature extraction step, the normal sound model updating step, the second acoustic feature extraction step, and the acoustic feature extractor updating step is repeatedly performed.
 20. An anomalous sound detection training method comprising: a first acoustic feature extraction step in which a first acoustic feature extraction unit extracts an acoustic feature of normal sound based on training data for normal sound by using an acoustic feature extractor; a normal sound model updating step in which a normal sound model updating unit updates a normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit; a second acoustic feature extraction step in which a second acoustic feature extraction unit extracts an acoustic feature of anomalous sound based on training data for anomalous sound and extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor; and an acoustic feature extractor updating step in which an acoustic feature extractor updating unit updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit, wherein processing at the first acoustic feature extraction step, the normal sound model updating step, the second acoustic feature extraction step, and the acoustic feature extractor updating step is repeatedly performed.
 21. A non-transitory computer-readable medium storing computer-readable instructions thereon, which, when executed by a computer, cause the computer to function as the units of the apparatus according to claim
 15. 22. An anomalous sound detection apparatus which detects whether or not input sound is anomalous sound, the anomalous sound detection apparatus comprising: an acoustic feature extraction unit that extracts an acoustic feature of the input sound based on an acoustic feature extractor; a degree-of-anomaly calculation unit that calculates a degree of anomaly from the acoustic feature of the input sound; and a determination unit that determines whether or not the input sound is anomalous sound based on the degree of anomaly and a threshold, wherein the acoustic feature extractor is optimized, by relating a probability density function followed by the normal sound to a probability density function followed by second anomalous sound, so that the acoustic feature extractor outputs an degree of anomaly to be determined as anomalous sound based on the threshold in case of normal sound and the acoustic feature extractor outputs an degree of anomaly to be determined as anomalous sound based on the threshold in case of the anomalous sound.
 23. The anomalous sound detection apparatus according to claim 22, wherein the probability density function followed by the normal sound is obtained from a first number of normal sound data, the probability distribution of the degree of anomaly followed by the second anomalous sound is obtained from a number of anomalous sound data sampled, the number of anomalous sound data sampled being substantially same as the first number of normal sound data.
 24. The anomalous sound sampling apparatus comprising, an anomalous sound sampling unit that generates an anomalous sound data where am amplitude spectrum of the anomalous sound data is based on normal sound, differs for each of anomalous sound data, is derived to be sparse and is smoothed in the temporal frequency direction, and an average spectrum of stationary noise is added to the amplitude spectrum of the anomalous sound data.
 25. The anomalous sound detection apparatus according to claim 23, wherein the anomalous sound data is generated by an anomalous sound sampling unit of an anomalous sound sampling apparatus comprising, an anomalous sound sampling unit that generates an anomalous sound data where am amplitude spectrum of the anomalous sound data is based on normal sound, differs for each of anomalous sound data, is derived to be sparse and is smoothed in the temporal frequency direction, and an average spectrum of stationary noise is added to the amplitude spectrum of the anomalous sound data. 